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Q\ • In recent years the presence of antiprotons (p's) in the cosmic ray (CR) flux incident upon the Earth has been 
firmly established by a series of balloon-borne experiments The CR p/p ratio has been measured from kinetic 

O ,■ energies of ~ 100 MeV to ~ 19 GeV. The measurements are summarized in Table 1. The CR p/p spectrum can be 
predicted using the Leaky Box Model (LBM) poT-pT2[ , which assumes that the p's originate from proton interactions in 
the interstellar (IS) medium. The p's then propagate within the Galaxy until they "leak out" with the characteristic 
CR Galactic storage time of about 10 million years (Myr) p3| . To obtain a prediction for the energy-dependent p/p 
ratio at the Earth, the IS spectrum must be corrected for modulation of the p and p fluxes as the particles propagate 
through the heliosphere. Thus, a comparison of the observed p/p spectrum with predictions tests the Galactic CR 
production-propagation model and the relevant heliospheric parameters. In addition, if the p's have a lifetime Tp- short 
compared to the Galactic storage time, the observed p spectrum will be depleted flp|Jl0|, particularly at low energies. 
This facilitates an interesting test of CPT invariance, which requires Tp — r p , where the proton lifetime t p is known 

" " ' to exceed 0(1O 32 ) yr @. 

In this paper the observed CR p/p spectrum is shown to be well-described by the predictions of an improved variant 
of the LBM, once modulation corrections are included. The data constrain the acceptable heliospheric properties. 
Extending the LBM to permit an unstable p, we obtain a limit on Tp significantly more stringent than current 
laboratory bounds obtained from searches for p decays in ion traps Jl~5|] and storage rings Q. 

In the LBM the only source of IS p's is the spallation of CR p's ]l0|-|l2|]: pNz — > pX, where Nz is a nucleus of 
charge Z, and X is anything. Realistic elemental abundances of the nuclei in the IS medium are given in ref s. JlT| , ^2| . 
The cross section for Z = 1 (the dominant contribution) has been measured in accelerator experiments ]10|. For 
Z > 1, some assumptions about nuclear physics are required to make use of the Z = 1 data. We have used the 
"wounded nuclcon" picture of Jl2| . The p's are assumed to propagate within the Galaxy until they are lost by one 
of two processes. The dominant loss process is leakage into intergalactic space. This process also applies to other 
CRs, including protons. The storage lifetime associated with this leakage mechanism varies with rigidity P and can 
be inferred from measurements of unstable secondary CR nuclei with known decay halnives. The sub-dominant p 
loss mechanism is pp annihilation. The annihilation rate is considerably smaller than the leakage rate. Leakage and 
annihilation, together with elastic and inelastic scattering of p's, are the ingredients of the standard leaky box model 
(SLBM), in which production and loss are assumed to be in statistical equilibrium. 
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Our analysis is based on the SLBM of Gaisser and Schaeffer using the parameters of E^] , but with a Galactic storage 
time improved to account for the non- uniform Galactic CR distribution [jl3| . This inhomogeneous LBM (ILBM) results 
in a somewhat longer storage time, 4.1 x 10 14 [1 + P/(3 GeV)] -0 6 sec. The uncertainties on the parameters of the 
ILBM result in uncertainties on the normalization of the predicted p/p ratio but, to a good approximation, do not 
introduce significant uncertainties in the shape of the predicted spectrum. Specifically, uncertainties on four ILBM 
parameters must be considered: (i) the storage time (±67% Q), (ii) the IS primary p flux (±35% Q), (hi) the p 
production cross section (±10% JlO|,p^[ ), and (iv) the composition of the IS medium, which introduces an uncertainty 
of < 6% on the predicted p flux Jl2|. In the following we neglect the last of these uncertainties since it is relatively 
small. Within the quoted fractional uncertainties on the other three parameters we treat all values as being a priori 
equally likely. Note that the predicted p/p ratio is approximately proportional to each of the three parameters under 
consideration. 

The curves in Fig. 1 show the ILBM p/p spectra for the parameter choices that result in the largest and smallest 
p/p predictions. The predicted IS spectrum does not give a good description of the observed spectrum at the top 
of the atmosphere. Good agreement is not expected because the CR spectra observed at the Earth are modulated 



as the particles propagate into the Sun's magnetic region of influence, the heliosphere 17 lq], which consists of the 



solar magnetic field B and the solar wind. Measurements by the Pioneer 10 and Voyager 1 spacecraft |19 indicate 
that the outer boundary of the heliosphere (the heliopause) is at heliocentric distance ro > 71 astronomical units 
(AU). Theoretical considerations favor ro = 100-130 AU jl7|]. Fortunately, the heliospheric corrections to the IS 
fluxes are not sensitive to variations of ro within these ranges. The wind, which is assumed to blow radially outwards, 
has an equatorial speed Vw ~ 400 km sec" 1 . Away from the equatorial plane the Ulysses spacecraft has found 
Vw ~ 750 km sec -1 Q. The wind pressure declines until it reaches the IS value at the heliopause. The solar wind 
plasma (nuclei and e~) moves in bulk outwards from the Sun with turbulence that peaks around solar magnetic 
maximum. The solar wind imposes an energy threshold for IS CRs reaching the Earth. This energy cutoff is much 
higher for particles coming in along the polar directions where the wind speeds are higher, than along the equatorial 
route where the wind speeds are lower. The wind also carries the solar magnetic flux outward. Solar rotation 17© twists 
the field lines to form a Parker spiral. The smoothed heliomagnetic field (£?© ~ 5 nT at the Earth's orbit) declines 
as it changes from radial at the Sun to azimuthal in the outer Solar System. The heliomagnetic polarity (sign(A)) is 
opposite in northern and southern solar hemispheres and switches sign somewhat after sunspot maximum (roughly 
every 11 years), when the field becomes more disordered. The regions of opposite magnetic polarity are separated by 
an approximately equatorial, unstable neutral current sheet. The sheet is wavy and spiraled; its waviness is measured 
by its "tilt" angle a, which relaxes from ~ 50° at polarity reversal to < 10° just before reversal p^7| , p0| -p2[ . 

Cosmic rays enter the heliosphere on ballistic trajectories and are then subject to the forces associated with the wind 
and solar mag netic field Hg. The propagation of the CRs within the heliosphere is described by a drift-diffusion 
(Fokker-Planck) equation [p3|. The CRs are pushed outwards by the bulk motion of the wind (elastic scattering), lose 
energy as they perform work on the wind (adiabatic deceleration or inelastic scattering), are diffusively scattered by 
field turbulence, and execute a drift orthogonal to the curving magnetic field lines as they spiral inwards along the 
field lines. Particles with qA > (< 0) drift in along a polar (sheet) route. The IS particles with sufficient energy to 
overcome the various energy losses reach the inner Solar System after being randomly scattered by turbulence. 

We compute the modulation of CR fluxes by the method of characteristics and combined Runge-Kutta/Richardson- 
Burlich-Stoer techniques |24|]. The calculation uses the heliospheric transport models of Jokipii et al. [S, combined 
with the ILBM p spectrum, and updated by Pioneer, Voyager, and Ulysses heliospheric measurements (l9|20[|. Our 
calculation includes magnetic curvature drift since older heliospheric models [^5| that neglected this drift component 
have been shown p2| to be inadequate. Note that magnetic drift changes sign for oppositely-charged particles and 
field polarities, in contrast to diffusion k which is charge-invariant and whose effect partly cancels in the p/p ratio. 
The calculation is simplified by ignoring turbulence where its effects are small, which in practice means everywhere 
except across the sheet, where for particles with speed v, we use: kj_ = (2 — 3) x 10 17 [£?©/-B(r)](P/GeV)°' 3 (u/c) 
m 2 /sec [|(|. Finally, IS "pickup" ions and solar CRs have been neglected [l7||20[| . 

In the following, we restrict our analysis to the data sets recorded by the MASS91, IMAX, BESS, and CAPRICE 
experiments (Table 1). These data were recorded in the period 1991-1994, corresponding to a well-behaved part of 
the solar cycle for which the heliospheric modulation corrections can be calculated with some confidence. In addition, 
compared to the earlier data, these more recent data sets benefited from detectors with improved particle identification 
and hence lower backgrounds. In Fig. 2 the measurements are compared with the predictions for a stable antiproton, 
taking the central parameter values for our ILBM calculation and three different sets of heliospheric corrections (lower, 
central, and upper). The corresponding heliospheric parameters are (a) Lower: equatorial Vw — 375 km sec" 1 , polar 
Vw = 700 km sec -1 , B§ — 4.0 nT; (b) Central: equatorial Vw = 400 km sec -1 , polar Vw = 750 km sec -1 , 
£>© = 4.5 nT; and (c) Upper: equatorial Vw — 425 km sec -1 , polar Vw = 800 km sec -1 , £?© = 5.0 nT. Note that 
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for each parameter set there are four curves, corresponding to the calculated heliospheric corrections at the dates 
of each of the four balloon flights. Given the variation of the predicted ratios with heliospheric parameter choice, 
we conclude that the predictions are able to give a good description of the data, and that there is no significant 
evidence for a short antiproton lifetime. Note that (i) for a given set of parameters, the heliospheric corrections are 
similar for the BESS, IMAX, CAPRICE, and MASS91 data, with appreciable differences only at the lowest kinetic 
energies; and (ii) in the region above ~ 3 GeV, the differences between the predictions corresponding to the different 
choices of heliospheric parameters are small compared to the statistical errors on the data. However, in the lowest 
energy region, below ~ 700 MeV, there are significant differences in the predictions for the three different choices of 
heliospheric parameters. Hence, in the fitting procedure employed to extract an upper limit on r-p, we must allow the 
overall normalization of the predictions to vary to account for the uncertainties arising from the choice of the ILBM 
parameters, and the heliospheric parameters to vary to account for the sensitivity of the low-energy p/p ratio to these 
parameters. 

To obtain a limit on Tp- we add to the ILBM one additional loss mechanism, p decay. Note that the p lifetime must 
be < O(10) Myr to significantly distort the p spectrum. However, this timescale is much longer than the lifetime 
sensitivity of existing laboratory searches for p decay Jlj|. The results from maximum likelihood fits to the BESS, 
IMAX, CAPRICE, and MASS91 measurements are shown in the inset of Fig. 2 as a function of the assumed Tp for eight 
heliospheric parameter sets Fl - F8. The Fl parameters are equatorial (polar) Vw — 375 (700) km sec -1 , and -B® = 
4.0 nT. The values for the equatorial (polar) Vw and -B® increase for the other parameter sets in steps of respectively 
5 (10) km sec -1 and 0.1 nT. Hence, the parameters for F8 are equatorial (polar) Vw = 410 (770) km sec -1 , and .B® 
= 4.7 nT. The fits, which take account of the Poisson statistical fluctuations on the number of observed events and 
the background subtraction for each data set (Table 1), also allow the normalization of the ILBM predictions to vary 
within the acceptable range (Fig. 1). The best maximum likelihood fits correspond to the heliospheric parameter sets 
F4 and F5, and are consistent with a completely stable antiproton. The extreme parameter sets (Fl and F8) are 
excluded at > 90% C.L. for all values of Tp-, and hence the data constrain the heliospheric parameters. In particular, 
the following parameter ranges are preferred: equatorial Vw = 375-410 km sec -1 , polar Vw = 700-770 km sec -1 , 
and £>® = 4.0-4.7 nT. Finally, the fits do not favor values of Tp- < 0(1) Myr. We obtain the bounds: 

Tp > 1.7 Myr (90% C.L.) , 1.1 Myr (95% C.L.) , 630 kyr (99% C.L.) . (1) 

A new generation of balloon experiments, or a dedicated orbital CR detector such as the recently tested AMS [ p7[ , 
could greatly improve the statistical precision of the measured p spectrum. This would facilitate more significant 
bounds on the heliospheric parameters, and enable a search for exotic sources of p's, particularly at low kinetic 
energies K. However we note that, neglecting diffusion, the solar wind cuts off the polar-routed flux at -ftT® = 0.26- 
0.40 GeV (mass-independent cutoff rigidity scaling approximately as gB® V w ). Particles with these threshold energies 
at the Earth have IS energies K\g — 0.54-0.80 GeV. In contrast, the sheet-routed flux exhibits no significant cutoff. 
Noting that in the present solar cycle the protons reach the Earth via the polar route, and p's reach the Earth via the 
sheet route, we conclude that when the next solar cycle begins (around 2003) new experiments will only be sensitive 
to IS p's with K m > 0.5 GeV. 

In summary, with an inhomogeneous LBM of IS p production and propagation, together with a calculation of 
heliospheric corrections that includes up-to-date data on the solar wind, we obtain good agreement with the observed 
CR p/p ratio. Extending the ILBM model to permit a finite Tp- we obtain lower limits on p decay that are more 
stringent than laboratory bounds. Our fits also constrain the ranges of the heliospheric parameters. We note that 
in the next solar cycle the effective cutoff in detectable low energy IS p's will be significantly higher than at present, 
which will decrease the sensitivity of experiments searching for exotic sources of antiprotons. 
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Experiment 


Field 
Pol. 


Flight 
Date 


KE Range 
(GeV) 


Cand- 
idates 


Back- 
ground 


Observed 
pip Ratio 


predict- 
tion 


Golden* 




+ 


June 1979 


5.6 - 12.5 


46 


18.3 


(5.2 ± 1.5) x 10~ 4 




Buffing'ton 1 " 


25 


n/a 


June 1980 


0.13 - 0.32 


17 


3.0 


(2.2 ±0.6) x 10~ 4 




Bogomolov 1 
Bogomolov* 
Bogomolov* 


-j 


_l_ 


1Q79 1Q77 
1984-1985 
1986-1988 


9 n - 5 n 
0.2 - 2.0 
2.0 - 5.0 


2 
1 
3 




(R 4- A\ v 1 n~ 4 
y\j m a ±\j 

(6 +14 ) x 10~ 5 

(2 4+?' 4 ~) x 10~ 4 
V^ ,4: — 1.3/ u 




PBAPJ 
PBAR f 
PBARt 


2£ 
3C 
3( 


- 


Aug. 1987 
Aug. 1987 
Alio- 1987 


0.205 - 0.64 
0.10 - 0.64 

n CiA _ 1 c;s 
yj .u^ i.jo 


n/a 
n/a 


n/a 
n/a 


< 4.6 x 10 _b 

< 2.8 x 10~ 5 

< 6.1 x 10 -5 


- 


LEAP 1 " 


31 




Aiip- 1987 


n 1 9 - n 86 


n/a - 


nTa - 


< 1.8 x 10~ 5 




MASS91 


§ 


+ 


Sep. 1991 


3.70-19.08 


11 


3.3 


(1.24^;-) x 10- 


1.3 x 10~ 4 


IMAX 
IMAX 
IMAX 


s 

s 


+ 
+ 
+ 


July 1992 
July 1992 
July 1992 


0.25 - 1.0 
1.0 - 2.6 
2.6 - 3.2 


3 

8 
5 


0.3 
1.9 
1.2 


(3.14+^t) x 10_b 
(5.36±|:|) x 10" 5 
(1.94+J'f) x 10" 4 


1.2 x 10 -5 
6.2 x 10~ 5 
1.1 x 10 -4 


BESS 


Z 


+ 


July 1993 


0.20 - 0.60 


6 


~ 1.4 


(5.2^-) x 10-" 


6.4 x 10~ b 


CAPRICE 
CAPRICE 


s 

9 


+ 

+ 


Aug. 1994 
Aug. 1994 


0.6 - 2.0 
2.0 - 3.2 


4 
5 


1.5 
1.3 


(2.5^) x W~ a 
(1.9+i'g) x 10- 4 


3.3 x 10 _b 
1.1 x 10~ 4 



TABLE I. Summary of cosmic ray antiproton results. Listed from left to right are the experiments, solar cycle polarity 
(note: minimum-to-minimum solar cycles 20:1964.6-1976.5, 21:1976.5-1986.7, and 22:1986.7-1996.3), balloon flight date, p 
kinetic energy range, number of p candidates observed, estimated number of background events, measured p/p ratio at the top 
of the atmosphere, and the ILBM prediction using the heliospheric parameters: equatorial (polar) Vw = 395 (740) km sec -1 , 
= 4.4 nT. f Not shown in Fig. 1 or used in analysis. J Not used in analysis. 
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FIG. 1. Observed p/p ratio at the top of Earth atmosphere (see Table 1). The solid curves show the upper and lower 
interstellar ratios predicted by the ILBM described in the text, without solar modulation. The broken curve shows the ILBM 
prediction with the same parameters used for the modulated predictions of Fig. 2. 
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FIG. 2. Observed p/p ratio at the top of Earth atmosphere compared with the ILBM predictions (see Fig. 1) after heliospheric 
modulation. The upper, central, and lower families of curves correspond to the upper, central, and lower heliospheric parameter 
sets described in the text. For each heliospheric parameter set there are four curves corresponding to the epochs of the four 
balloon flights (at low energies, from top to bottom: MASS91, IMAX, BESS, CAPRICE). The inset shows the fit results as a 
function of the assumed Tp for the eight heliospheric parameter sets (Fl - F8) described in the text. 
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